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Abstract. - We show that, in contrast with the usual electron-phonon Migdal-Eliashberg 
theory, the critical temperature T c of an isotropic s-wave nonadiabatic superconductor is strongly 
reduced by the presence of diluted non-magnetic impurities. Our results suggest that the 
recently observed T c -suppression driven by disorder in K3C60 (Phys. Rev. B 55, 3866 (1997)) 
and in Nd2-a ! Ce; C Cu04-,5 (Phys. Rev. B 58, 8800 (1998)) could be explained in terms of 
a nonadiabatic electron-phonon coupling. Moreover, we predict that the isotope effect on 
T c has an impurity dependence qualitatively different from the one expected for anisotropic 
superconductors. 



High-T c superconductors are narrow band systems with Fermi energies (Ep) one or two 
orders of magnitudes smaller than those of conventional superconductors (T) . In these materials 
therefore it appears unavoidable to address the validity of Migdal's theorem Q applied to the 
electron-phonon (el-ph) interaction and eventually other mediators 0. In fact, in fullerene 
compounds, bismuth oxides and high-T c cuprates, the typical phonon frequencies ujq can be 
comparable to Ep, making the quantity Xujq/Ef no longer negligible also for moderate values 
of the el-ph coupling constant A. Given this situation, the ordinary Migdal-Eliashberg (ME) 
theory, which is based on Migdal's theorem, hardly can appropriately describe the properties 
of such materials and a more general treatment of the problem should be employed. 

In general, by moving from the ME regime one could end to qualitatively different situations 
depending on the values of A and uio/Ep. For example, very strong el-ph couplings favour the 
formation of polarons and eventually of bi-polarons. This picture is certainly beyond the 
ME regime, however, as recently discussed in ref.j|, it is implausible that it can be at the 
basis of the phenomenon of high-T c superconductivity. Actually, one could go beyond ME 
regime without ending to polaron (bi-polaron) states by considering quasi-free charge carriers 
(A < 1) coupled nonadiabatically (0 < ujq/Ep < 1) to the lattice vibrations. In the following, to 
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distinguish this regime from the ME and the polaronic ones we use the concept of nonadiabatic 
fermions, which we define as quasiparticles (weakly) interacting nonadiabatically with phonons. 
In practice, such a regime can be formulated perturbatively by treating Xloq/Ef as the small 
parameter of the theory. At the zeroth order, the theory coincides with the ME limit while 
for finite values of Xloq/Ef the nonadiabatic fermions display anomalous behaviors, at the 
same time being faraway from the crossover to polarons. This latter feature is ensured by 
considering values of A smaller than the critical coupling A c (which is of order unity) of the 
crossover to the polaronic state. The reliability of such a perturbative approach is suggested 
by the comparison with exact results for the one electron system 0] and quantum Monte Carlo 
calculations for the many electrons case Q . 

The results obtained by previous studies of the nonadiabatic regime interest both the 
superconducting transition and the normal state properties. For example, the inclusion of 
the first vertex corrections beyond Migdal's limit have led to the possible enhancement of 
T c with respect to the adiabatic ME case jjj. Another interesting outcome is certainly the 
nonzero isotope effect of the effective mass of the nonadiabatic fermions || . The latter result 
provides a possible interpretation of recent measurements on cuprates j^]. 

In this paper, we present results on the non-magnetic impurity effects on the critical tem- 
perature T c and its isotope coefficient ar c for an isotropic s-wave nonadiabatic superconductor. 
We show that disorder strongly affect the nonadiabatic corrections and change their analytic 
properties substantially. As a result, T c can be strongly lowered with respect to the pure case, 
in contrast therefore with the adiabatic ME theory which, according to Anderson's theorem 
[ flo| , predicts an insensitivity of T c with respect to the presence of weak disorder for an isotropic 
s-wave superconductor. 

Our results can be of particular interest in light of the disorder dependence of T c recently ob- 
served in the fullerene compound K3C60 and in the electron doped cuprate Nd2- ;c Ce x Cu04_ l 5 
[ fi"ij Both compounds have a well-established s-wave symmetry of the gap, hence the 

disorder induced suppression of T c cannot be explained in terms of d-wave pairing Jl3| . It is 
however true that anisotropics could lead to qualitatively the same effect also if the gap is 
nodeless jl4|). Our results therefore provide an alternative interpretation based solely on the 
nonadiabatic regime of the el-ph interaction. 

To have a first idea of how disorder affects the critical temperature of a nonadiabatic 
superconductor, let us consider the first vertex correction in the normal state depicted in fig.l. 
In the figure, the solid and the wiggled lines are electron and phonon propagators, G and D, 
respectively, while the filled circles represent the el-ph matrix element <?(q) for momentum 
transfer q. From the usual rules for diagrams, the resulting vertex function P{k + q, k) is: 

P(k + q,k)=^2 g 2 (k - k')D(k ~ k')G(k' + q)G(k'). (I) 

k' 

Here, ^2 k and k are short notations for — ano - k, iw m , respectively. 

In the presence of impurities, the electron propagators entering eq. (^) are given by 
G _1 (fc) = iuj m — e(k) + iTuj m /\uj m \, where e(k) is the electronic dispersion and T = l/2r 
is the impurity scattering rate [15| . We neglect for the moment the self-energy due to the el-ph 
coupling, the effects of the finite bandwidth on the impurity contribution and the momentum 
dependence of g(k — k'). By employing some approximations which will be specified later, the 
vertex function (Q) can be evaluated as a function of the dimensionless momentum transfer 
Q = q/2kp, where kp is the Fermi momentum, and the exchanged frequency u>. In fig. 2 
we show the sign of the vertex function in the Q-lo plane for the half-filling case and for 
ujq/Ef — 0.5. The different lines denote the boundaries where the vertex changes sign. On 
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the right side of the lines the vertex is positive while on the left side is negative. The effect of 



impurities is twofold. First, the non-analyticity in uj — 0, Q = found for the pure case |16 
(solid line) is removed when r ^ 0. Second, by increasing the value of T the boundary lines 
shift toward higher values of the exchanged frequency, reducing therefore the region where the 
vertex is positive. When we average the vertex over the exchanged momentum and frequency, 
we find that by increasing T the average becomes negative. Since the vertex correction (|lj) 
enters into the generalized Eliashberg equation for a nonadiabatic superconductor [jl6), we 
expect from the above result that T c should be lowered with respect to the pure case. 

In order to confirm this hypothesis, we solve the Eliashberg equations beyond Migdal's 
limit by including the effect of disorder in the Born approximation. The relevant diagrams 
for both the normal, Ejv, and anomalous, Eg, self-energies are depicted in fig. 3. In terms of 
the renormalized frequency W n — uj n Z n , where Z n = 1 — T,N(iuj n )/iuj n , and the renormalized 
s-wave gap function <p n — Z(iLu n )A(iuj n ), the generalized Eliashberg equations reduce to: 

4>n^n = irT c ^{\[l + 2\P{Q c] n,m)]D(uj n -Lo m ) + \ 2 C(Q c ;n,m)}-^r ]m , (2) 

c ; n, m)]D(u! n - uj m ) m rj m , (3) 

m ' 1 

where r) m = (2/ir) arctan(£ , i?/|W / m |) is the finite band- width factor and = 1 — ^nT/lWnl is 
the renormalization induced by the elastic impurities jnl. 

In writing the above equations, we have assumed a dispersionless phonon spectrum with 
frequency c^o so that the phonon propagator is given by: D(wi) = ^/[(iuji) 2 — ojq]. The 
terms P(Q c ;n,m) and C(Q c ;n 7 m) are the vertex and cross corrections, respectively |t], [l6|| . 
The parameter Q c = q c /2kp is an upper-cutoff over the momentum transfer which follows 
from the model we use for the el-ph coupling function g 2 (Q) — (g 2 / 'Q 2 C )8(Q C — Q) where 9 
is the Heaviside step function ]?], This model simulates the effect of the strong electron 
correlations which lead to mainly forward el-ph scattering processes [^). In fact, in a strongly 
correlated system, a single charge carrier is surrounded by the corresponding correlation hole 
whose linear size £ c can extend over many lattice units [H . The charge carrier can therefore 
respond only to charge modulations with wave vector q c < l/£ c - 

The solution of equations (d)-^) without the vertex and cross corrections P(Q C ', n, m) and 
C(Q c ;n,m) has been already reported by Choi in rcf . JTsf] . The main result was that, for 
finite values of Ep, T c can be slightly lowered by the presence of impurities. Here, we solve 
instead the completely self-consistent equations with the inclusion of both the vertex and 
cross corrections. In evaluating P(Q c ;n,m) and C(Q c ;n,m) we use a density of states Np 
constant over the entire electron bandwidth 2Ep and we employ a small momentum transfer 
approximation which is valid for small values of the dimensionless cutoff parameter Q c . From 
the above discussion, the small Q c approximation should be suitable for strongly forward el-ph 
couplings resulting from the effect of strong electron correlation. After the integration over 
the energy and the s-wave average over the momentum transfer Q have been performed, the 
vertex and cross corrections reduce to the following form: 



P(Q c ;n,m) = T c D(u) n - lji) j B(n, m, I) + 



2 




AEpQ 2 c 



A(n, m, I) - B(n, m, l)(Wi - Wi- n+m ) 
(2EpQl) 2 

02 \ 2 \ 1 

,(4) 
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C{Q c ;n, m) 



D(ui n — u>i)D(uji — uj m ) < 2B(n, —m, I) + arctan 

A(n, -to, I) - B(n, -to, l)(W t - Wi- n -m) 2 
X 2E F Q%\Wi — Wi- n - m \ 



4E F Q 2 C 



\Wi - Wi- n . 



(5) 



A(n,m,l) = (Wi-Wi- n+m ) 



arctan 



(w) 



arctan 



+m 



(6) 



B{n,m,l) = (Wi- Wi- n+m ) 



E F Wi 



— n+m 



E F 
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W? 
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(7) 



To obtain the critical temperature T c , we solve numerically the self-consistent equations 
(||)~(@) f° r different values of the adiabatic parameter ojq/Ef, the momentum transfer cut-off 
Q c and the impurity scattering rate T. The dependence of T c upon the adiabatic parameter 
is shown in fig. 4 for pure (T/uio — 0, solid lines) and disordered (T/u>o = 0.5, dashed 
lines) nonadiabatic superconductors for different values of Q c . Our self-consistent calculations 
confirm the results of ref. |?], i.e., small values of Q c lead to an enhancement of T c with 
respect to the Migdal limit. Moreover, at ujq/Ef = 0, T c is independent of V in agreement 
with Anderson's theorem ]Tc| |. On the other hand, when ujq/Ef > the impurities lower T c 
for all values of Q c and ujo/Ef- We interpret this behavior in terms of the enhanced negative 
contribution of the vertex and cross corrections induced by the presence of the impurities as 
depicted in fig. 2. This negative contribution leads to a reduction of the effective nonadiabatic 
el-ph pairing interaction resulting in the reduction of T c . 

In fig. 5 we show T c as a function of T for luq/Ef = 0.2 (a) and ujq/Ef — 0.4 (b). The 
thin solid lines refer to the case without vertex and cross corrections and correspond to the 
approximation scheme used in ref. fl9| . As seen also by the inserts of fig. 5, when we include 
the vertex and cross corrections (thick lines), the reduction of T c with the increase of T can 
be much stronger than the reduction given by only the finite bandwidth effects. 

So far, we have investigated the effects of disorder on a nonadiabatic superconductor with an 
s-wave symmetry of the order parameter. However, anisotropies of the gap lead to an impurity 
dependence of T c which, to a first approximation, can be described by using the Abrikosov- 
Gorkov (AG) scaling law p0[ modified in order to represent d-wave, anisotropic s-wave and 
other types of symmetries of the order parameter 13, 14J: ln(T c o/T c ) = x[ 1 i r (l/2+7)— ^(1/2)], 
where ^ is the digamma function and 7 = T/(2irT c ). The parameter x is a measure of the 
anisotropy of the order parameter: x = 1 (x — 0) for a d-wave (s-wave) superconductor Q. 
According to the AG law, for x 7^ the impurities induce a monotonous reduction of T c in a 
way qualitatively similar to the one observed for the nonadiabatic case. 

Given the above situation it could be therefore difficult to decide whether the T c suppression 
observed in K3C60 Q and in Nd 2 _ a; Ce a ;Cu04_5 (l2| should be ascribed to anisotropies of the 
order parameter or instead to the nonadiabatic el-ph interaction. Here, we propose that a more 
suitable quantity to look at could be the ion-mass dependence of the critical temperature. In 
fact, the isotope coefficient ax c resulting from the AG-type relation is: 



"T C 



d\n(T c /T< 
G?ln7 



c:0 ) 



(8) 
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where aT a0 is the isotope coefficient for the pure system. Since T c /T c o decreases by increasing 
7, equation (||) predicts a monotonous impurity induced enhancement of a.T a compared with 
the corresponding value in the pure limit ]2lf| . Such a behavior is qualitatively different from 
the one displayed by an isotropic s-wave nonadiabatic superconductor. In fact, as it is shown 
in fig. 6, where we report numerical results for for the same parameters as in fig. 5(a), the 
decrease of T c /T c0 due to the impurities is accompanied by a non-monotonous dependence of 
(Xt c at least for small values of the momentum cut-off Q c . Moreover, for larger values of Q c 
the isotope coefficient decreases with the impurity concentration showing therefore a behavior 
opposite to the one given by eq.@. It would be therefore important to measure the isotope 
coefficient and its evolution with the amount of disorder in the two s-wave superconductors 
K3C60 and Nd2- 2: Ce :I ;Cu04_5. Such a measurement could in fact decide whether the observed 
T c suppression in these materials JlT| |l2| is given by anisotropy or by the nonadiabatic regime 
of the el-ph interaction. We note moreover that a measurement of ar c vs T c /T c o could provide 
an experimental tool for an estimation of the typical momentum scattering Q c in addition to 
the one obtained by tunneling measurements [ p2[ . 

In summary, we have shown that in the nonadiabatic regime the critical temperature is 
lowered by non-magnetic impurities which also lead to an unusual impurity dependence of 
the isotope coefficient. Our results together with recent measurements on K3C60 and 
Nd2-xCe x Cu04_5 |0| suggest that the el-ph interaction in these systems could be in the 
nonadiabatic regime. We propose also that in order to confirm or disregard this hypothesis, a 
suitable experiment could be the measurement of the isotope coefficient as a function of the 
amount of disorder. We conclude by noticing that, to our knowledge, there are no experimental 
results on the effect of disorder on the bismuth oxides. Since these materials are s-wave 
superconductors p3j with estimated loq/Ef — 0.14 p3], we predict that in these material T c 
should be considerably lowered by non-magnetic impurities and eventually display an isotope 
effect with anomalies as described above. 

C. G. acknowledges the support of a INFM PRA project (PRA-HTCS). 
REFERENCES 

[1] Uemura Y. J. et al: , Phys. Rev. Lett, 66 (1991) 2665. 

[2] Migdal A. B., Sov. Phys. JETP, 7 (1958) 996. 

[3] Schrieffer J. R., J. Low Temp. Phys., 99 (1995) 377. 

[4] Chakraverty B. K., Ranninger J., and Feinberg D., Phys. Rev. Lett., 81 (1998) 433. 

[5] Capone M., Ciuchi S., and Grimaldi C, Europhys. Lett., 42 (1998) 523. 

[6] Freericks J. K., Zlatic V., Chung W., and Jarrel M., Phys. Rev. B, 58 (1998) 11613. 

[7] Grimaldi C, Pietronero L. and Strassler S., Phys. Rev. Lett., 75 (1995) 1158. 

[8] Grimaldi C, Cappelluti E., and Pietronero L., Europhys. Lett., 42 (1998) 667. 

[9] Zhao G. M., Hunt M. B., Keller H. and Muller K. A., Nature, 385 (1997) 236. 
[10] Anderson P. W., J. Phys. Chem. Solid, 11 (1959) 26. 
[11] Watson S. K. et al: ., Phys. Rev. B, 55 (1997) 3866. 
[12] Woods S. I. et al: , Phys. Rev. B, 58 (1998) 8800. 
[13] Radtke R. J. et al: , Phys. Rev. B, 48 (1993) 653. 
[14] Openov L. A., JETP Lett, 66 (1997) 661. 

[15] Rickayzen G., Green's Functions and Condensed Matter, (Academic, New York) 1980. 
[16] Pietronero L., Strassler S. and Grimaldi C, Phys. Rev. B, 52 (1995) 1995; Grimaldi C, 
Pietronero L. and Strassler S., Phys. Rev. B, 52 (1995) 10530. 



6 



EUROPHYSICS LETTERS 



[17] Zeyher R. and Kulic M., Phys. Rev. B, 54 (1996) 8985; Grilli M. and Castellani C, Phys. 
Rev. B, 50 (1994) 16880; Keller J., Leal C. E., and Forsthofer F., Physica C, 206-207 
(1995) 739. 



[18] Danylenko O. V. et at: , cond- mat/9710234, Preprint 1997 



[19] Choi H. Y., Phys. Rev. B, 53 (1996) 8591. 

[20] Abrikosov A. A. and L. P. Gor'kov L. P., Sov. Phys. JETP, 12 (1961) 1243. 

[21] Bill A., Kresin V. Z., and Wolf S. A., Z. Phys. B, 104 (1997) 759. 

[22] Ummarino G. A. and Gonnelli R. S., Phys. Rev. B, 56 (1997) R14 279. 

[23] Brawner D. A., Mamser C, and Ott H. R., Phys. Rev. B, 55 (1997) 2788. 



M. SCATTONI et al: ANOMALOUS IMPURITY EFFECTS ETC. 



7 
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Fig. 1. - The electron-phonon scattering process. The filled circles represent g(q). The last 
diagram is the first vertex correction g(q)P(fc + q, k) which in the adiabatic limit gives a 
negligible contribution according to Migdal's theorem. 
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Fig. 2. - Sign of the vertex function for different values of the impurity scattering rate T at 
uq/Ef — 0.5. Solid line: T = 0; dashed line: T = 0.1u; ; dot-dashed line: T = 0.5o;o- 
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Fig. 3. - (a): electronic self-energy including the first vertex correction beyond Migdal's limit 
and the impurity contribution, (b): Self-consistent equation for the anomalous self-energy 
generalized to include the first nonadiabatic diagrams (vertex and cross) and the impurity 
contribution. 
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Fig. 4. - Critical temperature as a function of the adiabatic parameter u> /Ep for A = 0.7. The 
solid lines refer to the pure case (r = 0) while the dashed lines are the results for T — O.buo . 
Note that at w^/Ep = the critical temperature is independent of T. 
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Fig. 5. - Critical temperature as a function of the impurity scattering rate T for A = 0.7, 
u/Ep = 0.2 (a) and u>o/Ep = 0.4 (b). The thick (thin) lines refer to the case with (without) 
the nonadiabatic corrections (Q) and (||). Thick solid: Q c = 0.1, thick dashed: Q c = 0.3, 
thick dot-dashed: Q c = 0.5. Inserts: same calculations rescaled with respect to the critical 
temperature T c q of pure systems. 
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Fig. 6. - Isotope coefficient ax c as a function of T c in the presence of non-magnetic impurities. 
Both quantities are normalized to their corresponding values ar c0 and T c o for the pure limit. 
The curves refer to the case of fig. 5(a). 



